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1.  Summary  of  Research  Carried  out  in  this  Contract 

In  this  research,  the  goal  was  to  develop  techniques  on  how  to  enhance  the  received 
signals  in  a  near  field  multi-input  multi-output  (MIMO)  environment  where  beam  forming  is  not 
possible.  This  is  done  through  the  use  of  adaptivity  on  transmit.  This  technique  is  based  on  the 
principle  of  reciprocity,  is  independent  of  the  material  medium  in  which  it  is  transmitting,  and 
incorporates  near-field  environments  and  multipath.  The  objective  here  is  to  select  a  set  of 
weights  adapted  to  each  receiver  to  be  applied  to  each  transmitting  antenna,  which  is  a  function 
of  the  user  location,  so  that  the  transmitted  signal  at  the  carrier  frequency  may  be  directed  to  a 
particular  receiver  location  while  simultaneously  minimizing  the  received  signal  strengths  at 
other  receiver  locations. 

Next  the  principles  of  polarization  was  utilized.  The  polarization  adaptivity  on  transmit 
can  be  used  to  enhance  the  received  signals  directed  to  a  pre-selected  receiver  in  a  near-field 
multi-input  multi-output  (MIMO)  environment.  The  objective  here  is  to  select  a  set  of  weights  on 
the  transmitting  antennas  adapted  to  each  receiver  based  on  the  principles  of  reciprocity.  Using 
the  polarization  properties,  when  the  number  of  receiving  antennas  is  greater  than  the  number  of 
transmitting  antennas,  the  transmitted  signal  may  be  directed  more  to  a  particular  receiver 
location  while  simultaneously  minimizing  the  reception  signal  strength  at  other  receivers. 
Wideband  performance  of  this  system  has  also  been  studied. 

In  addition,  the  effect  of  errors  on  the  spatial  location  of  antennas  and  how  they  impact 
the  performance  in  adaptive  processing  for  antenna  arrays  including  its  effect  on  STAP  had  also 
been  studied.  This  final  report  deals  with  this  topic. 

2.  Introduction 

The  direct  data  domain  least  squares  (D  LS)  approaches  are  based  on  the  work  of  Sarkar 
and  Sangruji  [1].  Then  Sarkar  et  al  [3,  4]  presented  two  more  variations  to  this  approach,  the 
Backward  processor  and  the  Forward-Backward  processor.  These  algorithms  were  then  extended 
to  the  space-time  domain,  where  Park  [5]  described  a  generalized  eigenvalue  based  processor 
and  Sarkar  et.  al.  [6]  implemented  the  D  LS  processors. 
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The  direct  data  domain  least  squares  space-time  adaptive  processing  (D  LS  STAP) 
algorithm  has  been  used  for  suppressing  highly  dynamic  clutter  and  interference  [6-8].  This 
enables  the  system  to  detect  potentially  weak  target  returns.  In  this  analysis,  we  assume  that  the 
system  consists  of  a  linear  array  of  N  equally  spaced  antenna  elements.  We  assume  that  the 
system  processes  a  number  of  pulses  (M)  within  a  coherent  processing  interval  (CPI).  Each  pulse 
repetition  interval  (PRI)  consists  of  the  transmission  of  a  pulsed  wavefonn  of  finite  bandwidth 
and  the  reception  of  the  reflected  energy  by  the  radar  aperture  and  receivers,  with  a  bandwidth 
matched  to  that  of  the  radar  pulse.  Data  can  be  arranged  into  a  three-dimensional  matrix 
commonly  referred  to  as  a  data  cube  [7].  We  assume  that  the  signals  entering  the  array  are 
narrowband,  with  a  planar  phase  front,  and  consist  of  a  signal-of-interest  (SOI)  along  with 
interference  plus  noise.  The  noise  (thermal  noise)  originates  in  the  receiver  and  is  understood  to 
be  independent  across  elements  and  pulses.  The  interference  is  external  to  the  receiver  and 
consists  of  clutter,  jammers,  and  signal  multipath,  which  may  be  coherent  or  incoherent.  For  a 
uniformly  spaced  linear  array  (ULA),  the  complex  envelope  of  the  received  SOI  with  unity 
amplitude,  for  the  mth  pulse  and  nth  antenna  element,  can  be  described  by 

Sm n  =  exp[/2n{^  cos(^)  +  m  ^fd }],  for  m  =  1,2,  ...,  M;  n  =  1,2,  ...,  N  (1) 

k  lr 

where  d  is  the  inter-element  spacing  between  the  antenna  elements,  k  is  the  wavelength,  fr  is  the 
pulse  repetition  frequency,  </>s  is  the  angle  of  arrival  and  Ij  is  the  Doppler  frequency  of  the  SOI. 
For  the  nth  antenna  element  and  mth  pulse,  the  complex  envelope  of  the  received  signal  is  given 
by  the  scalar  quantity 

XM=  as  •  S  +  Interference  +  Noise  (2) 

where  as  is  the  complex  amplitude  of  the  SOI  entering  the  array.  Between  the  scalars  Xm  n  and 
Xm  n+i.  there  is  a  phase  difference  only  for  the  SOI,  which  is  given  by  [6,  8] 

Z,  =exp[/27r^cos(4)]  (3) 

A 

Also,  there  is  a  phase  difference  for  the  SOI  in  the  scalars  Xm+i,nand  Xm  n..  due  to  the  Doppler  of 
the  SOI,  which  is  given  by  [6,  8] 
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(4) 


Z2  =  exp(_/27iy-) 

For  the  two-dimensional  case  these  differencing  are  performed  with  elements  offset  in  space  due 
to  the  subscript  m,  and  in  time  due  to  the  subscript  n,  and  jointly  in  space  and  time  due  to  the 
subscripts  m+1,  n+1  [6,  8].  The  three  types  of  difference  equations  are  then  given  by 


Xm,n  -  z. 

•  V 

yvm,n+l 

(5) 

Xm,n  -  z;1 

•  V 

yvm+l,n 

(6) 

m.n  -  Z(‘ 

•  z2‘  •  xm+1>n+1 

(V) 

Now,  we  can  form  the  cancellation  rows  in  matrix  [F],  whose  elements  are  generated  using 
equations  (5)-(7)  [6,  8],  which  will  deal  with  only  the  undesired  signals,  as  the  SOI  has  been 
cancelled  out  of  equations  (5)-(7).  And  the  elements  of  the  first  row  of  matrix  [F]  is  implemented 
by  fixing  the  gain  of  the  array  along  the  direction  of  the  SOI  in  both  space  and  time  as 


[1,  Zj,  Zf,---,  ZjNa_1 ,  Z2 ,  ZtZ2 ,  Zf Z2 , •  •  •  Z^a_1Z2 ,  Zj ,  Z^-.Z^Z”'-1]  (8) 


where  Na  is  the  spatial  dimension  and  Nt  is  the  temporal  dimension  of  the  window.  Na  and  Nt 
must  satisfy  the  following  equations 

N+1 


N  < 


N,  < 


2 

M+1 


(9) 

(10) 


The  resulting  matrix  equation  is  then  given  by 

[F][W]  =  [Y]  =  [C  0  ...  0]T.  (11) 

where  [W]  is  the  two-dimensional  STAP  weight  vectors  and  C  is  the  scalar  gain  of  the  array 
along  the  SOI  in  space  and  time. 

The  second  process  is  to  use  Backward  processor,  which  can  be  implemented  by 
conjugating  the  element-pulse  data  and  processing  this  data  in  reverse.  We  also  use  Forward- 
Backward  process  that  is  formed  by  combining  the  Forward  and  Backward  solution  procedure  [3, 

4]- 
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When  a  D3LS  STAP  is  used,  we  assume  that  each  antenna  element  is  designated  to  be 
placed  at  a  specific  location  [8],  In  the  real  environment,  however,  it  is  possible  for  the  location 
of  the  antenna  elements  to  fluctuate  due  to  an  installation  error  or  environmental  effects. 
Research  on  random  perturbations  in  the  excitation  coefficients  of  the  antenna  arrays  has  been 
the  subject  of  several  papers  [9,  10].  These  papers  consider  random  perturbations  of  the 
amplitude  and  phase  of  the  amplifiers.  As  a  follow-on  to  the  paper  by  Zaghloul  [9]  on  the 
deterioration  of  the  equivalent  isotropically  radiated  power  (EIRP)  due  to  random  perturbations 
in  the  excitation  coefficients  of  the  antenna  array,  Hwang  and  Sarkar  [12]  present  a  paper  for 
allowable  tolerances  in  the  position  of  antenna  elements  for  acceptable  performance.  In  this 
paper  we  identify  how  random  position  errors  can  affect  a  D  LS  approach  when  applied  to 
STAP.  Our  goal  is  to  find  how  the  system  will  perfonn  when  the  antenna  elements  are  randomly 
moved. 

In  Section  II,  we  discuss  EIRP  degradation  of  the  array  antennas  due  to  random  position 
errors  in  terms  of  the  normalized  EIRP  degradation.  In  Section  III,  to  explore  the  utility  of  the 
EIRP  degradation,  a  probabilistic  analysis  of  a  simple  phased  array  antenna  is  studied.  In  Section 
IV,  the  D  LS  approaches  for  STAP  are  simulated  using  the  relationship  between  the  EIRP 
degradation  and  position  errors  due  to  random  locations  of  the  elements  in  an  array. 

3.  EIRP  Degradation  of  Array  Antennas  Due  to  Random  Position 
Errors 

Consider  an  antenna  array  consisting  of  N  antenna  elements  fed  by  N  separate  amplifiers. 
Except  for  a  scaling  factor,  the  far  field  E{6,  (/))  is  given  by 

E(0,<!>)  =  X  ^'(“x'+v^wz,) Eel{0,  </>)  (12) 

i=\ 

where 

a,  =  |  a,  |  ei¥i  complex  excitation  coefficient  for  the  /lh  antenna  element 
x,  =  the  random  perturbation  of  the  x-coordinate  of  the  zlh  antenna  element 
y;  =  the  random  perturbation  of  the  y-coordinates  of  the  zth  antenna  element 
z;  =  the  random  perturbation  of  the  z-coordinates  of  the  /th  antenna  element 
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u  =  ksm6cos(j),  v=  ksinOsin</>,  w=  k  cos#,  k=2n/X 
Eel(9,<j) )  =  the  antenna  element  pattern  along  the  direction  (6,(f>) . 

The  probability  density  functions  of  the  random  perturbations  x„  v;,  and  z,  may  be 
uniform  or  have  Gaussian  distributions.  For  uniform  distributions,  the  probability  density  of  x„ 
Vi,  and  z,  are  given  by 


P(x,)  = 


=  0, 


xi-Ax<xi<xi+Ax 

otherwise 


(13) 


p(yt)  = 


yl-Ay<yl<yl+Ax 


=  0,  otherwise 


(14) 


P(z,)  =  ,  z,~  Az  <  z,  <  zi  +  A  v 

=  0,  otherwise 

For  Gaussian  distributions,  the  probability  density  functions  for  x„  v„  and  z,  are  given  by 

/  \  ^  -(Xj—x,)2/ 2ctT- 

P(x1)=  , —  e 


yflnt 


na,. 


p(y,)  = 


na. 


/  x  1  -(z  -z,-)2  /2cr2 

P(Zi)  =~i=  g 


V27 t< 


na_ 


(15) 


(16) 

(17) 

(18) 


where  x, ,  v,  and  z,.  are  the  mean  values  of  the  x,. ,  y, ,  and  z,  -positions  respectively,  and  a],  a2y 


and  a2  are  the  variance  of  the  x,  y  and  z-coordinates  respectively. 

According  to  Zaghloul  [9],  EIRP  along  the  direction  (&,</>)  is  proportional  to  the  power 
radiated  in  the  same  direction.  Therefore  EIRP,  is  given  by 


J  =  1 


N  N 

i= 1  7  =  1 


g.g: 


in  which  *  denotes  the  complex  conjugate,  and 


(19) 
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Gj  =  at  ■  Eel  {6,  <!).  (20) 

The  probability  distribution  of  the  EIRP  will  be  Gaussian  with  the  mean  value  of  P 
denoted  by  P  and  a  standard  deviation  of  <j(P)  from  the  central  limit  theory  [11].  The 
probability  distribution  of  the  EIRP  is  within  the  confidence  interval 
[p-na(P)]/P0,  [P  +  n<j{P)\l P^\  .  Left  side  of  the  confidence  interval  in  the  probability 

distribution  of  the  EIRP,  [P  -na(P)\/  P0 ,  represents  the  worst  case  for  all  possible  cases. 

Normalizing  P  and  <j(P)  with  respect  to  the  unperturbed  power  P0  presents  the  degradation  in 

the  EIRP  of  the  array.  This  quantity,  \P  -  ncr(P)]/P0 ,  is  the  normalized  EIRP  degradation. 

Along  the  same  line  of  the  amplitude  and  phase  perturbation  in  reference  [9],  we  choose  n  =  3  so 
that  the  worst-case  degradation  of  the  EIRP  is  within  99.7  percent  of  all  possible  cases  for  the 
given  random  position  distributions. 

The  expected  value  P  and  the  standard  deviation  <j(P)  of  P  can  be  expressed  in 
separable  terms  containing  the  statistical  parameters  for  jq,  v,  and  z,-.  Each  amplifier  is  located  at 
each  x,  y,  and  z-position  and  is  considered  to  be  independent.  They  have  the  same  distribution. 
An  expression  for  the  influence  of  random  locations  of  the  antenna  elements  can  be  calculated  by 
integrating  the  product  of  the  probability  distribution  and  the  power  P.  The  expected  value  P  is 
obtained  as 

N  N  N 

P  =  YJ\Gif+  Hx  Hv  H.  Y  Y  ejMx,~"j)+vif,~fj)+w{Zi~Zi)]  Gi  G*  (2 1 ) 

i=i  i=i  j= i 

where  Hx  is  the  x-position  variation  dependent  factor,  Hv  is  the  y-position  variation  dependent 

factor  and  H:  is  the  z-position  variation  dependent  factor.  These  expressions  are  given  in  Table 
1. 
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Table  1:  x,  y  and  z-position  Factors  in  the  Expected  Value  Expression. 


Factor 

Uniform  Distribution 

Gaussian  Distribution 

Hx 

(sinwAv/wAv)2 

2  2 
—U  <JY 

e 

Hy 

(sinvAv  /  vAy)2 

2  2 
-V  (7 

e  y 

K 

(sinwA_  /  wA_)2 

2  2 
-W  <7t 

e 

The  next  step  in  calculating  the  standard  deviation  <j{P)  indicates  the  extent  of  the  array 
EIRP  degradation.  The  standard  deviation  is  given  by 

a(P)  =  yJvav(P)  (22) 

and  the  variance  is  defined  as 


var(P)  =  E[P2]-P2 

in  which  E[-]  denotes  the  expectation  value.  We  obtain  EfP2]  as 


(23) 


N  N 


N  N 


E[/>!]  =  Z|G,|  +2x£I|G,|  Gj  +4 xH,,HIyH,,XZe 


j[u(xi -Xj  )+v(yi -yj  )+w(zi -z}  )]  \q\2QQ* 


i= 1 


i= 1  y=l 


N  N  N 


i=l  y=l 


+  4  x  HlxHlyHXz 


J\M  (*i  -x,  )+v(y,  -y,  )+w(z,  -z, )] 


G,\  GG 


l  ] 


+  H3xHlyHXzYJ 


*=1  1=1  7=1 
k^i^j 

N  N 

j2[u(Xj  -Xj  )+v(yf  -yj  )+w(zf  -zy. )]  ^2  ^*2 

/=1  7=1 
i*j 

N  N  N 


j[u(xk  +xt -2 Xj  )+v(yk  +yt -2y j  )+ w(zk  +zt -2 zy  )]  ^  ^  ^*2 


£=1  i=l  7=1 
k^j 


_l_  ^7 tM (“**  +2*,-  j  )+v(—yk  +2^/ ->> y  )+ w(-z*  +2z,- -z y  )]  J 


N  N  N  N 


k= 1  /=1  1=1  7=1 


j[u  (xk  -xt  +Xf  -Xj  )+v(yk  -y,  +yt  -yj  )+w(zk  -z,  +zf  -z  - )] 


*  I  *  ./ 


(24) 


where  the  factors  //, ,  ,  II,  ,  and  II, .  result  from  estimating  constituents  of  the  second  moment 
of  the  power  P,  and  are  given  in  Table  2  for  the  Uniform  and  Gaussian  distributions. 
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Table  2:  x,  y  and  z-position  Factors  in  the  Variance  Expression. 


Factor 

Uniform  Distribution 

Gaussian  Distribution 

Hi,x 

sinwA^ 

V  UA*  J 

2 

2  2 
— tt  cr„ 

£ 

Hl,y 

f  •  .  \ 

sinvAr 

vAv 

V  y  J 

2 

2  2 
-V  <7 

Hi,z 

sin  wAz  'l 

l  ^Az  J 

2 

-wV 

e 

H2,X 

/  N 

sin  2nA  v 
2u  A 

V  X  y 

N2 

^t«2a2 

e 

H2,y 

Z  N 

sin  2vA  v 

2vA 

X  y  y 

-4v2ct2 

e  y 

h2jZ 

/ 

V 

sin2wAz 

2wA. 

y 

/ 

-4w2a? 

e 

H3jx 

^sinwA^ 

2 

f  sin  2«A  v  N 

-3  w2ct2 

e 

V  “A.v  j 

l  2“Ax  J 

H  3jy 

sinvA,,  'l 

y  , 

2 

f  sin  2vAv  ^ 

y 

-3v2ct; 

e 

vA  „ 

V  y  ) 

i  2vA-  j 

h3jZ 

sin  wA z  | 

2 

f  sin2wAz 

-3w2<7^ 

e 

.  wA--  J 

^  2wAz 

H4jX 

sinuA^ 

l  wA.  J 

4 

-2  u2a; 

e 

H4,y 

sinvA,, N 

l  vA.v  J 

4 

T  2  2 

-2v  <jv 

e  y 

H4,z 

sin  wAz^l 

l  ^Az  J 

4 

-2w2  <j2 

e 
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4.  Example  of  EIRP  Degradation  In  Antenna  Arrays 

To  explore  the  utility  of  the  EIRP  degradation,  a  probabilistic  analysis  on  a  simple  linear 
phased  array  antenna  as  shown  in  Figure  1  is  performed.  Using  this  example,  we  investigate  how 
the  EIRP  degradation  behaves  as  a  function  of  the  various  parameters. 


O 


n=,l 


d=0.4775  A. 


y 

- o  - x 

n=10 

At  z  =  0  plane 


Figure.  1:  The  Uniformly  spaced  linear  array  configuration 

Consider  a  linear  array  consisting  of  10  elements  as  shown  in  Figure  1 .  The  array  consists 
of  10  isotropic  omni-directional  point  radiators  and  is  fed  by  10  amplifiers.  The  elements  of  the 
linear  array  are  equally  spaced  and  their  inter-element  spacing  along  the  x-direction  is  given  by 
0.4775A  (a=1  m),  where  X  denotes  the  wavelength  of  operation.  To  simplify  the  example,  we 
assume  that  all  the  elements  are  located  at  z  =  0  so  that  the  z-positions  of  all  the  elements  are 
fixed.  Therefore  the  EIRP  degradation  of  the  linear  array  is  caused  by  random  perturbations  in 
the  x  and  y  position  which  can  have  two  types  of  probability  distribution,  either  uniform  or 
Gaussian.  Hz  is  1  in  the  Expected  value  of  eq.(21)  and  Httz  is  also  1  in  the  second  moment  of 

the  power  of  eq.(24).  When  the  positions  of  all  the  elements  are  unperturbed,  Hx  and  Hv  equal 

unity,  and  also  P  becomes  the  unperturbed  power  P0 .  As  we  expect,  the  variance  is  identically 
zero,  as  given  by  eq.(23). 

The  expected  value  and  the  standard  deviation  now  is  applied  to  the  linear  array  case  to 
get  the  normalized  EIRP  degradation  with  n  =  3  and  \P  -  3 cr(E’)]/  P0 .  Of  special  interest  is  the 
EIRP  degradation  within  the  half  power  beam  width  between  (j)  =  90°  and  <j)  =  95.33°.  Figures  2 
and  3  show  the  worst  case  degradation  with  3o(P).  Figure  2  represents  the  case  for  the  Gaussian 
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perturbation  with  ax  and  y  =  0.0,  0.03,  0.05,  and  0.07  [m],  and  Figure  3  deals  with  uniform 
perturbation  with  A  x  and  y  =  0.0,  0.03,  0.05,  and  0.07  [m].  As  expected,  the  unperturbed  case 
shows  a  0  dB  value  for  all  the  angles.  The  EIRP  degrades  with  the  larger  standard  deviation  of 
the  Gaussian  perturbation  and  the  uniform  bounds  for  a  uniform  perturbation.  The  EIRP  due  to 
the  random  positions  of  the  antennas  is  significantly  decreased. 


Linear  Array:  Gaussian  Distribution 


Figure.  2:  Normalized  EIRP  degradation  with  n  =  3  for  Gaussian  perturbation 
with  standard  deviation  of  axandy  =  0.0,  0.03,  0.05,  and  0.07  [m] 

The  normalized  EIRP  degradation,  [P  -3cr(P)]/  P0,  represents  the  worst  case  for  99.7 
percent  of  all  possible  cases  as  mentioned  before.  This  value,  however,  depends  on  different 
values  of  a  for  a  Gaussian  distribution  or  has  a  different  A  for  a  uniform  distribution,  as  we  see 
from  Figures  2  and  3.  The  question  is  what  is  the  allowable  tolerance  for  a  and  A  and  for  the 
elements  to  operate  appropriately  in  the  antenna  system  for  adaptive  signal  processing.  We  focus 
on  the  main  lobe,  especially  inside  the  half  power  beam  width.  To  choose  the  appropriate  axandy 
and  Axandy,  we  select  the  point  that  makes  the  3dB  loss  below  the  optimum  for  the  normalized 
EIRP  degradation  graph,  at  the  angle  corresponding  to  the  half  power  beam  width.  At  this  -3  dB 


10 


point,  the  antenna  performance  using  conventional  beam  fonning  will  be  significantly  decreased. 
However,  the  3dB  loss  of  the  normalized  EIRP  degradation  at  the  angle  of  the  half  power  beam 
is  the  allowable  worst  case  for  an  antenna  system  to  operate  properly  compared  to  the 
unperturbed  case.  From  these  a  and  A  obtained  from  the  3dB  loss  point  of  the  normalized  EIRP 
degradation,  we  investigate  how  the  system  behaves  if  the  antenna  elements  are  randomly 
perturbed  by  those  a  and  A. 


Figure.  3:  Normalized  EIRP  degradation  with  n  =  3  for  uniform  perturbation 
with  different  bounds  of  Axandy  =  0.0,  0.03,  0.05,  and  0.07  [m]. 

As  illustrated  in  Figure  4,  the  -3dB  point  of  the  normalized  EIRP  degradation  with  n=3 
of  a  uniform  linear  array  occurs  at  95.33°  with  axandy  =  0.042  m  and  Axandy  =  0.072  m.  Figure  4 
also  illustrates  that  the  relationship  between  the  normalized  EIRP  degradation  with  n=3  and  ctx  and 
y,  and  between  the  normalized  EIRP  degradation  and  A  x  and  y  at  a  specific  angle,  95.33°.  The 
larger  the  axandy  and  Axandy,  the  lower  the  normalized  EIRP  degradation.  We  apply  this  value, 
gx  and  y  =  0.042  m  and  A  x  and  y  =  0.072  m,  and  then  investigate  how  the  random  position  errors 
affect  the  D3LS  approach  for  STAP. 
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Linear  Array,  at  angle  95.33 


Figure.  4:  The  relationship  between  the  normalized  EIRP  degradation  with 
n  =  3 ,  and  ax  and  y  and  A  x  and  y  at  the  angle,  95 .33°  for  a  linear  Array. 

5.  Simulation  Results 

The  received  signals  modeled  in  this  simulation  consist  of  the  SOI,  main  beam  clutter, 
discrete  interferers,  jammers  and  thermal  noise.  The  clutter  is  modeled  as  point  scatterers  placed 
approximately  every  0.1  degrees  apart.  The  amplitude  of  the  clutter  points  is  modeled  with  a 
normal  distribution  about  a  mean  that  result  in  a  signal  to  clutter  ratio,  SCR,  of -12.7  dB.  Ten 
strong  point  scatterers  are  modeled  in  this  simulation,  based  on  the  angle  of  arrival  (AO A)  and 
Doppler  parameters  as  defined  in  Table  3.  Thermal  noise  generated  in  each  receive  channel  is 
independent  from  one  another.  The  resulting  SNR  is  approximately  30  dB.  The  jammer  is 
modeled  as  a  broadband  noise  signal  that  arrives  from  100°  in  azimuth,  and  covers  all  Doppler 
frequencies  of  interest.  Summing  the  power  of  the  interfering  sources  received  by  the  first 
channel  and  comparing  it  to  the  power  of  the  SOI,  the  average  input  signal  to  interference  plus 
noise  ratio  (SINR)  is  evaluated  as  -18.6  dB.  And  random  position  errors  of  the  antenna  element 
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are  axandy  =  0.042  m  and  Axandy  =  0.072  m.  While  processing  the  signal,  the  antenna  elements 
are  randomly  moved.  In  this  simulation  we  consider  two  situations.  One  in  which  the  antenna 
elements  at  every  pulse  repletion  interval  (PRI)  sampling  instance  have  different  spatial  positions 
from  the  previous  time  instance.  Another  situation  is  that  the  antenna  random  locations  are  fixed 
within  a  coherent  processing  interval  (CPI). 


Table  3:  Parameters  Related  to  the  Simulation. 


Wavelength 

1  m 

Signal 

AOA 

0 

O 

Pulse  Rep.  Freq 

4  kHz 

Doppler 

1169  Hz 

Number  of 

Antennas 

N=  10 

SNR 

30  dB 

Number  of  Pulses 

M=  16 

Jammer 

AOA 

100° 

Half  Power 

Beam  width 

10.66° 

Doppler 

Covers  all  Doppler 

frequencies  of  interest 

Process¬ 

or 

Forward 

Na  =  7,  Nt  =  9 

Discrete 

interferes 

AOA 

[  85°  120°  40° 

35°  30°  140° 

100°  70°  50°  125°] 

Forward- 

Backward 

Na  =  8,  Nt  =  9 

Clutter 

extent 

Main  beam 

Doppler 

[400  -800  1700 

-1400  325  -1650 

950  -1200  -125  1450]  Hz 

Point 

scatterers 

0.1°  apart  in 

angle 

3 

With  N  =  10  channels  and  M  =  16  pulses,  the  performance  of  the  three  D  LS  algorithms, 
namely  the  Forward,  Backward  and  the  Forward  and  Backward  method  is  evaluated  based  on  the 
number  of  adaptive  weights  and  the  output  SINR.  The  Forward  and  the  Backward  methods 
utilize  7  spatial  (Na)  and  9  temporal  (Nt)  degrees  of  freedom  (DOF)  resulting  in  a  total  of  63 
DOF,  while  the  Forward-Backward  method  employs  8  spatial  and  9  temporal  DOF,  for  a  total  of 
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72  DOF.  The  spectrum  of  the  input  signal  is  shown  in  Figure  5.  Here  the  circles  indicate  the 
locations  of  the  discrete  interferers  and  the  triangles  define  the  location  of  the  main  beam  clutter 
and  the  +  mark  indicates  the  location  of  the  SOI. 


2000 
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-500 

-1000 

-1500 

-2000 


0  20  40  60  80  100  120  140  160  180 


Angle  (Deg) 


Figure.  5:  The  spectrum  of  the  Input  Signal. 

Table  4  shows  the  output  signal-to-interference-plus-noise-ratio  (SINR)  for  the  unperturbed 
case,  and  perturbed  case  with  a  Gaussian  profile,  and  a  uniform  profile.  These  values  are 
averaged  over  100  runs.  The  resulting  weights  for  the  unperturbed  case  and  for  the  perturbed 
case  with  Gaussian  and  uniform  density  functions  using  the  Forward  method  are  shown  in  Figure 
6,  7  and  8,  when  the  antenna  elements  at  every  time  instance  have  different  spatial  positions  from 
the  previous  time  instance.  Antenna  positions  are  randomly  located  with  axandy=  0.042  mor  Ax 
and  y  =  0.072  m.  The  Output  SINR  of  a  Gaussian  and  uniformly  perturbed  cases  are  lower  than  the 
Output  SINR  of  the  unperturbed  case.  Even  though  we  choose  different  a  and  A,  we  get 
approximately  an  output  SINR  of  ±1  ldB.  This  shows  that  we  get  EIRP  degradation  which  can  be 
obtained  analytically.  From  Figures  7  and  8,  the  system  with  random  positions  generates  nulls 
slightly  moved  from  the  angle  of  arrival  (AOA)  of  discrete  interferers  as  compared  to  the 
unperturbed  case.  The  output  SINR  of  the  unperturbed  case  is  about  18-19  dB  and  the  output 
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SINRs  for  the  Gaussian  and  uniform  density  distributions,  which  are  about  1 1  dB  when  the 
antenna  elements  at  every  time  instance  have  different  spatial  positions  from  the  previous  time 
instance.  This  perturbed  system  still  removes  the  jammer,  clutters,  discrete  interferers  and  noise. 
When  the  random  locations  of  the  antenna  are  fixed  in  a  CPI,  one  also  obtains  almost  the  same 
output  SINR.  The  value  of  the  output  is  dependent  on  the  value  of  the  tolerance  of  the  antenna 


Table  4:  The  output  signal  to  interference  plus  noise  ratio. 


Case  A 

(When  antenna  random  locations 
are  changed  at  every  time  instance) 

Case  B 

(When  antenna  random 
locations  are  fixed  in  a  CPI) 

Forward 

Unperturbed 

18.210  [dB] 

Gaussian 

11.458  [dB] 

11.028  [dB] 

Uniform 

11.480  [dB] 

10.515  [dB] 

Backward 

Unperturbed 

18.210  [dB] 

Gaussian 

11.458  [dB] 

11.028  [dB] 

Uniform 

11.480  [dB] 

10.515  [dB] 

Forward- 

Backward 

Unperturbed 

19.213  [dB] 

Gaussian 

11.481  [dB] 

11.150  [dB] 

Uniform 

11.844  [dB] 

11.183  [dB] 

elements. 

From  these  results,  as  long  the  antenna  element  position  errors  are  within  the  limits  set  by  a 
=  0.042  m  for  the  Gaussian  distribution  and  A  =  0.072  m  for  the  uniform  distribution,  the 
adaptive  array  has  a  reasonable  perfonnance  for  STAP.  If  we  want  to  provide  stricter  tolerances 
in  the  distribution  of  the  antenna  elements,  i.e.,  a  <  0.042  m  for  the  Gaussian  distribution  and  A 
<  0.072  m  for  the  uniform  distribution,  then  one  can  get  higher  output  SINR  which  may  be  close 
to  the  output  SINR  for  the  unperturbed  case.  As  a  result,  one  can  predict  how  much  the  system 
would  be  degraded  by  analytically  solving  the  problem  using  the  nonnalized  EIRP  degradation. 
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Unperturbed  Case  by  Forward  Method  jg 


0  20  40  60  80  100  120  140  160  180 

Angle  (Deg) 

Figure.  6:  Weight  Spectrum  of  the  Forward  Method  for  the  unperturbed  case. 


Gaussian  Case  by  Forward  Method  dB 


0  20  40  60  80  100  120  140  160  180 

Angle  (Deg) 

Figure.  7:  Weight  Spectrum  of  the  Forward  Method  for  the  perturbed  Gaussian  profile 

with  CJX  and  y  =  0.042  m. 
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Uniform  Case  by  Forward  Method 
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Figure.  8:  Weight  Spectrum  of  the  Forward  Method  for  the 
perturbed  Uniform  profde  with  Axandy  =  0.072  m. 

6.  Conclusion 

In  this  report,  expressions  are  derived  for  the  expected  value,  standard  deviation,  and  the 
nonnalized  EIRP  degradation  with  n=3  due  to  the  random  locations  of  the  array  elements.  The 
nonnalized  EIRP  degradation  represents  the  worst  case  for  99.7  percentile  of  all  possible  cases. 
We  investigate  the  effects  of  random  antenna  positions  on  a  D3LS  approach  for  STAP  using 
these  EIRP  degradations  according  to  the  random  position  of  the  array  elements.  Even  though  a 
=  0.042  m  for  the  Gaussian  distribution  and  A  =  0.072  m  for  the  unifonn  distribution  are  given 
for  the  antenna  elements,  we  can  still  get  an  acceptable  output  SINR.  When  the  antenna  elements 
at  every  time  instance  have  different  spatial  positions  from  the  previous  time  instance,  the  output 
SINR  is  almost  the  same  as  in  the  case  where  the  antenna  elements  are  located  at  fixed  random 
locations  within  a  CPI.  In  this  report,  the  results  for  the  EIRP  degradation  have  been  predicted 
for  the  Gaussian  and  for  the  uniformly  perturbed  cases.  Random  antenna  position  errors  degrade 
the  output  SINR  of  a  D3LS  Approach  for  STAP. 
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